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A theory of the collapse of a punctured antibubble is developed. The motion of the rim of 
air formed at the edge of the collapsing air film cannot be described by a potential flow and is 
characterized by high Reynolds numbers. The rim velocity is not constant but gradually decreases 
with time and is determined by the balance between the surface tension and hydrodynamic drag 
forces. A collapse equation is derived and solved. The agreement between the theory and existing 
experiments is shown. 
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An antibubble is a thin spherical gas film containing 
and being surrounded by a liquid. It is a peculiar an¬ 
tipode to an ordinary bubble, a spherical liquid film with 
a gas inside and outside. The study of antibubbles has a 
relation to the physics of films, interfaces, foams, bubbles, 
and drops [l[. The formation of antibubbles differs from 
the cavitation of bubbles: antibubbles are generated by 
gently dripping or pouring a surfactant solution onto the 
surface of the same solution. Antibubbles were seemingly 
observed in 1931 during the investigation of soap drops 
on a water surface 0: some of the drops sank under 
the water surface and showed interference colors; from 
this observation the authors concluded that they had ob¬ 
served drops surrounded by a soap film. The existence 
of the air film was mentioned later Q, and the objects 


studied were referred to as “inverted soap bubbles” [JJ or 
“inverse bubbles” . The term “antibubble” was coined 
by Pavlov-Verevkin in his Russian paper entitled “Soap 
Antibubbles” in 1966 Interest in antibubbles was 
revived in 2003, when the formation and collapse of an¬ 
tibubbles and the development of the concomitant fluid 
instabilities were observed with a high-speed video cam¬ 
era 0 . Since then the stabilization of antibubbles d, 0 , 
optimal conditions for the antibubble formation [l^ , and 
antibubble lifetime distribution [lll - [l^ have been stud¬ 
ied. The formation of an “antidrop,” an object similar to 
an antibubble but in which a liquid phase replaces the air 
film, was also observed 15|- Antibubbles not only allow 
one to observe and study new fast microhydrodynamic 
phenomena but also are interesting by their potential ap¬ 
plications EH- 


Recently, a new experimental study of the collapse of 
an antibubble was conducted [l8j |. in which an attempt 
was made to find the factors that determine the veloc¬ 
ity of the edge of the shrinking air film after puncturing 
the antibubble with a pin. The authors concluded that 
this velocity is virtually constant. However, the earlier 
experimental data Q definitely show that the velocity 
decreases during the collapse. This surprising paradox 
poses a question about our understanding of the antibub¬ 


ble collapse and requires a proper theoretical description 
of the phenomenon, still absent despite the existing ex¬ 
periments. The purpose of this Letter is to develop the 
theory of the antibubble collapse. In particular, the the¬ 
ory solves the above paradox, reveals what actually de¬ 
termines the collapse velocity, and shows that the an- 
tibubble collapse differs from the rupture of liquid films. 

Figure [T] shows a scheme of a collapsing antibubble. 
The antibubble of thickness d and radius R is p unctured 
at point P at time t = 0. In the experiment [l8| the char¬ 
acteristic values for the antibubble thickness and radius 
are d ~ 3 ^m and i? Ri 4 mm, respectively; therefore, we 
will assume that d <C R. After puncturing, the liquids 
inside and outside the antibubble come into contact with 
each other at point P and the air film, originally closed, 
starts shrinking due to the surface tension of the liquid. 
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The film shrinking results in the appearance and expan¬ 
sion of a circular hole in the antibubble. The collapsing 
air film is axisymmetric with respect to polar axis OP, 
where O is the center of the antibubble, so that the cen¬ 
ter O' of the hole always remains at OP while moving 
from starting point P to diametrically opposite point Q. 

The edge of the hole represents a rim that contains all 
the air from the collapsed part of the antibubble. We 
assume, for simplicity, that the cross section of the rim 
by the plane passing through OP is a circle of radius a; 
then the rim is a torus. Consideration of the rim makes 
sense ii a ^ d. We also assume that a <C i?, which can 
be rewritten as ttmax ^ R, where ttmax = y2(iH?/n is 
the maximum rim radius that corresponds to the torus 
of maximum volume 2Tr'^a^^y. = AirR^d. 

The position of the rim is given by polar angle 0, the 
angle between polar axis OP and radial vector OP', where 
P' is a point at the guiding circle of the torus. Equating 
the volume of the rim, V-dm = 27r^a^i?sin0, to the volume 
of the collapsed air film, Woi = 27ri?^d(l — cos6>), yields 
the dependence of the rim radius on 0: 


a = aoYtan-, (1) 

where oq = \JdRl'K is the semicollapse rim radius, the 
rim radius corresponding to 0 = 7''/2, when the first 
hemisphere of the antibubble has disappeared. The vol¬ 
ume conservation condition is utilized in Refs. 0, m 
to estimate the thickness of the antibubble; the pos¬ 
sibility of using this condition is experimentally shown 
in Ref. [I^. Equation o makes sense for 0 from the 
range 0min ^ 0 < 0max, where 0min = ^d/irR and 
(^max = TT — tt^^x/R- The first inequality is necessary 
for the possibility of neglecting the difference between 
using 0 and the more correct value 0 + a/R in calculat¬ 
ing VcoV, it is in fact equivalent to the aforementioned 
inequality d. The angle 0max in the second inequal¬ 
ity corresponds to the torus of maximum volume, so that 

a(0 max) — Q-max- 

The rim moves in the liquid with velocity v = Rd9/dt. 
We assume that the rim is sufficiently thin: a <C Rs'md. 
If this condition is satisfied, the flow about the rim is 
quasi-two-dimensional. This condition implies the above 
condition a R and can be violated only at the final 
stage of the collapse, when the polar angle is close to tt. 
Thus, we should assume that tt — 0 ^ a/R, or equiva¬ 
lently TT — 0 TT — 0niax; this means that 0 should not be 
very close to 0max- 

Denote by p the mass density and by rj the dynamic 
viscosity of the liquid. We can then define the corre¬ 
sponding Reynolds number Re = 2avp/r]. Experimen¬ 
tally [l^, p ~ 10^ kg m“^, rj ^ 10“^ Pa s, w ^ 1 ms“^, 
and a can be estimated via d and i? as a ^ oq 60 /rm; 
therefore. Re ^ 120. The flow about the cylinder at such 
Re is characterized by the existence of the von Karman 


vortex street [l9| and is obviously nonpotential due to 
the existence of vortices. 

Note in this connection that the assumption about the 
flow potentiality was used in Ref. 18| and allowed the 
authors, after some order-of-magnitude estimations, to 
obtain a constant velocity of the rim, v oc yja/pd, where 
cr is the surface tension. We conclude that the potential 
flow about the rim is not realized and that the above 
formula is inapplicable in the experimental conditions of 
Ref. 0. Particularly, v need not be constant, which re¬ 
solves the paradox mentioned in the introduction. Thus, 
the antibubble collapse is characterized by a nonpoten¬ 
tial, vortical flow with high Reynolds numbers, and it is 
necessary to study which factors actually determine the 
velocity of the rim. 

The motion of the rim is due to the surface tension. 
The surface tension force acting upon the unit segment 
of the rim is = 2 (t, with 2 due to the existence 
of two interfaces between the air film and liquid. This 
force is counterbalanced by the hydrodynamic drag force 
Fd = Cdpv^a, where Cd is the drag coefficient. For sim¬ 
plicity we will assume the drag coefficient constant, which 
implies relatively high Reynolds numbers; we may choose 
Cd ~ 1.1—the typical dra g co efficient for a circular cylin¬ 
der at 300 ^ Re < 2 X 10^ [20l|. We also restrict ourselves 
to considering only the two mentioned forces, and the 
comparison with experiments will show a posteriori that 
other possible forces do not make a significant contribu¬ 
tion. 

The above balance between the forces can be quali¬ 
tatively interpreted as thus: the surface energy of the 
collapsing air film turns into the kinetic energy of the 
liquid behind the rim. If the unit segment of the rim has 
moved the unit distance, the energy of the disappeared 
air film is 2cr whereas the volume of the liquid is 2a; the 
kinetic energy in this volume is ^ pv^a because the typ¬ 
ical energy density is ^ pv^/2. The balance between the 
two energies is equivalent to the condition Fo- = Fd with 
Cd ~ 1. Such consideration makes sense only for Re 1. 

From Fcr = Fd with use of Eq. CD we obtain the rim 
velocity as a function of 0: 


V = Vo 



where 


I'D 


2g 

Cdpao 


( 2 ) 

( 3 ) 


is the semicollapse rim velocity, which corresponds to 
0 = 7r/2. We immediately see from Eq. ([2|) that the rim 
velocity is not constant and decreases with time because 
0 obviously increases during the collapse. 

Let us discuss an interesting question: Why is the 
rim velocity found in Ref. [l8| constant if the hydro- 
dynamic drag force vanishes in a potential flow due to 
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d’Alembert’s paradox while the surface tension force is 
nonzero? The authors consider the balance between the 
rate of increasing the kinetic energy of the potential flow 
about the rim and the rate of decreasing the surface en¬ 
ergy. Since the energy of potential flow about the cylin¬ 
der of unit length is maU^/2, where rua = p'KO^ is the 
added mass (2]|, the former rate results in the action 
of an added mass force. This force is, however, not the 
usual added mass force TOaW due to acceleration, which 
vanishes, but the reactive force ihaV due to increasing 
the added mass of the rim. That is why the result of 
Ref. [l^ , V oc ajpd^ is similar to that of Ref. : in 
such a consideration the air film behaves similarly to a 
liquid film with a massive rim. 


Actually, the situation is quite different: the flow 
is nonpotential, and the hydrodynamic drag force is 
not only nonzero but also much larger than the added 
mass force. Both the velocity and the added mass 
change with time, so we estimate the latter force as 
Fa = dma,v/dt. It follows from Eqs. ([T]) and ([2]) that 
Fa/Fd = (37r/4C'd)(a/Fsind) <C 1. Therefore, the col¬ 
lapse of an air film is not similar to the collapse of a liquid 
film. 


Let us define a typical time characterizing the collapse, 
T = R/vq, and a dimensionless time r = t/T. From 
Eq. m we then get the collapse equation 


do J e 

* = V“‘2’ 


(4) 


where 9 is considered as a function of r. Thus, Eq. 
determines the time dependence of 9. 

First find the semicollapse time to, the time of the 
disappearance of the first hemisphere of the antibub¬ 
ble. Integrating •^tan(d/2) over 9 from 0 to 7r/4 yields 
the dimensionless semicollapse time tq = /3(5/8) « 
1.209, where j3{x) = [tpHx + l)/2) — 'ip{x/2)]/2, ijj^x) = 
d In r(a;) / dx is the digamma function, and r(a;) is the Eu¬ 
ler gamma function [23j. The semicollapse time is then 
to = tqF. 


Second find the collapse time tcoi, the time of the 
disappearance of the whole antibubble. Integrating 
•ytan(0/2) over 9 from 0 to tt/2 yields the dimension¬ 
less collapse time TcoI = 7rsec(7r/8) ~ 3.400 [i^. The 
collapse time is then tcoi = TcoiT. We see again that the 
collapse gradually slows down: tcoi/to ss 2.812 > 2. 

Now we turn our attention to the time dependence 
of 9. With the substitution y{9) = ■^tan(d/2) we derive 
an analytical solution of the collapse equation [Eq. (jT])] 



X 


FIG. 2. (color online) Polar angle 6 against dimensionless 
time r—exact solution [Eq. ([5])] of the collapse equation 
[Eq. ([4])]. Inset: absolute (a.e.) and relative (r.e.) errors 
of fit [Eq. ([Til] against r. 


in the form r = t(9): 


r(0) = cos — ^ < 2 arctan 


n=0 


y{9) CSC — -I- (—1)" cot — 


+(-!)” In 


y{9f + {-l)-2y{9)sm- + l 


— Sin ■ 


n—0 


< 2 arctan 


y{9) sec — -I- (—1)" tan - 


+(-l)’^ln 


y(0)2 + (_l)n2y(0)cos- + l 


• (5) 


This solution is presented in Fig. [5] 

The exact solution [Eq. ([5])] is rather cumbersome, and 
we need a tractable analytical fit to easily work with the 
inverse function 9 = 9{t). We see directly from the col¬ 
lapse equation [Eq. (|4])] that 9 has asymptotics 


p(5T/8)4/5, r^O, 

6 »(t) = <1 7r/2-b (r - To), t^tq, (6) 

[tT - 2 [3(Tco1 - t)/ 8]'‘^^, T^Tcol. 


We then construct a fit 9{t) to the exact solution 9{t) 
from the condition that 9{t) satisfies the exact asymp¬ 
totics [Eq. (|6])]: 

rAoT‘‘/®(l - Ar“), O^r^To, 

9{t) = W “ Fo(tco1 - t)"*/^ (7) 

[ X[1 -b B{Tco\ - t)^], 

where Aq = 2(5/8)^/^ 1.373, = 2(3/8)'^/^ 

0.5408, a = (27r/5 - to){Aot^^^ - tt/2)-^ ^ 1.703, 
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FIG. 3. (color online) Polar angle 6 (A) and angle velocity 
dO/dt (B) against time t —experimental data (balls and 
circles) and theoretical curves (solid lines). 


P = (tcoI - To - 27r/3) [7r/2 - So(tco 1 - tq)^/^] ^ « 3.069, 
A = (1 - 7r/2AoT^/®)r-“ r; 1.260 x IQ-^, and B = 
[7r/2Bo(Tcoi - To)'‘/^ - l](rcoi - to)“^ Ri 1.846 x 10“^. 
The inset in Fig. [2] shows the accuracy of the fit: the 
absolute error 0 — 0 < 4 x 10““^ ~ 0.02° and the relative 
error [6 — 6)/9 < 0.02%. We see that the constructed 
fit [Eq. Q] is of high accuracy and hence can be used 
instead of the exact solution [Eq. ([5])]. 

Figure [3] shows the data from Ref. Q on the time 
dependence of polar angle 9 and the theoretical curve 
9{t/T) (A). We observe a good agreement between the 
theory and experiment. Least squares fitting gives the 
characteristic time T = 18.06 ± 0.08 ms (standard error 
shown), from which we calculate the semicollapse time 
to = 21.8 ± 0.1 ms and collapse time tcoi = 61.4 ± 0.3 ms. 

Figure [3] also shows the data from Ref. [3] on the time 
dependence of angle velocity d9/dt and the theoretical 
curve {cot[0(t/r)/2]}^/^/r (B). We again observe a good 
agreement, and we can independently calculate T from 
these data: T = 16.8 ± 0.5 ms. This value is consis¬ 
tent with the above value, but has a much higher uncer¬ 
tainty because the experimental error of d9/dt is much 
higher than that of 9. These errors can be estimated 
from the respective fits: erg ~ 0.02 and (Jde/dt ~ 7 
and correspond to the typical relative errors Se ~ 1.3% 
and 5de/dt ~ 12% (calculated at the semicollapse point, 
where 9 = n 12 and d9/dt = T~^). Thus, due to smaller 
experimental errors, the data on the time dependence of 9 
contain more precise information about various temporal 
characteristics of the collapse. 

Figure |4] shows the data from Ref. [l^ on the time 
dependence of rim velocity v and the theoretical curve 
uo{cot[0(t/T)/2]}^/^, where vq = R/T and R = 3.80 mm 
(A). From least squares fitting we have T = 4.1 ± 0.2 ms 


FIG. 4. (color online) Rim velocity v (A) and distance be¬ 
tween antibubble and hole centers, r, over radius R (B) 
against time t —experimental data [1^ (triangles and balls) 
and theoretical curves (solid lines). Inset: r and R on the 
image [l^ for the third data point. 


and the semicollapse rim velocity vo = 0.93 ±0.05 ms“t^ 
and we estimate the absolute and relative experimental 
errors as ~ 0.2 ms“t and Jt, ~ 21%. These data are 
consistent with the fact that the rim velocity gradually 
decreases with time and do not require that it is constant. 

Reference [1^ also contains a chronological sequence 
of images of the collapsing antibubble [Fig. 2 therein]. 
I process these images in a manner shown in the inset 
in Fig. m to obtain the data on the time dependence of 
r/i?, where r is the distance between the center O of the 
antibubble and the center O' of the circular hole in the 
air film (Fig. |T]). Note that the rim is not clearly visible 
in the images because its small radius is much less than 
the antibubble radius. These data, with error bars from 
irregularity of the air film edge, and the theoretical curve 
cos9{t/T) are presented in Fig. S] curve (B), and again 
demonstrate a good agreement between the theory and 
experiment. Weighted least squares fitting gives more 
precise values than above: T = 3.9 ±0.1 ms and vq = 
0.97 ±0.03 ms-b 

Very significantly, we can independently find the semi¬ 
collapse rim velocity vq directly from Eq. ©: using the 
experimental values d = 2.95 fxm, p = 992 kgm“^, and 
a = 0.033 Nm“^ 18| and putting Cd = I.l, we obtain 


Vq = l.OI ± 0.05 ms“^, which is consistent with the two 
above values. This consistency confirms that the tempo¬ 
ral development of the collapse observed in the experi¬ 
ment is correctly described by the theory, without fitting 
but from the known physical parameters. 

When choosing Cd, we have considered the drag force 
acting upon a solid cylinder, with no-slip boundary con¬ 
ditions. Meanwhile, the rim is made of air, and free-slip 
boundary conditions may be more suitable, but result in 
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remarkably decreasing Cd [2^. On the other hand, we 
have the air rim not in pure water but in a soap solution, 
and surfactant molecules adsorb at the interface. Because 
of the film of surfactant molecules, an ordinary bubble 
experimentally behaves at Re <C 1 not as an air ball with 
free-slip conditions but as a solid ball with no-slip condi¬ 


tions [2^, l26j , and a similar effect may be expected in the 
case of the rim. However, the effect can be reduced be¬ 
cause the boundary of the moving rim goes into the two 
fixed interfaces of the air film before the rim. Thus, it is 
reasonable to consider Cd as an adjustable parameter. 


We have seen that Cd 1 agrees with the experimen¬ 
tal data and hence could argue in favor of no-slip condi¬ 
tions. Interestingly, the theory allows us to calculate Cd 
from experiment with the help of Eq. ®. I propose to 
study experimentally how Cd depends on various phys¬ 
ical parameters of the antibubble and on the surfactant 
concentration, which can shed some light on the problem 
of boundary conditions. 


In conclusion, I have developed the theory of the col¬ 
lapse of an antibubble. After the antibubble is punc¬ 
tured, the air film starts shrinking due to the surface 
tension of the liquid, which results in the appearance 
and expansion of a circular hole in the antibubble. The 
collapse is characterized by a nonpotential flow about the 
rim that forms at the edge of the shrinking air film, and 
the flow has high Reynolds numbers. The rim velocity 
does not vary linearly with yajpd and is not constant; 
therefore, the collapse of an antibubble differs from the 
rupture of a liquid film 22|. The velocity gradually de¬ 
creases with time and is determined by the balance be¬ 
tween the surface tension and hydrodynamic drag forces 
acting upon the rim. The collapse occurs so that the sur¬ 
face energy of the collapsing air film turns into the kinetic 
energy of the liquid behind the rim. I have derived and 
solved the collapse equation, the solution of which gives 
the time dependence of the polar angle and describes the 
temporal development of the collapse. I have demon¬ 
strated the agreement between the theory and currently 
available experimental data. In regard to an interesting 
problem of single bubble sonoluminescence , I pro¬ 

pose to study how the antibubble collapse occurs in the 
presence of an acoustic field. 
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